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The role of negative buoyancy in convective Cepheid 
models. Double-mode pulsations revisited 
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CN ; Abstract 

C"| \ The longstanding problem of modeling double-mode behaviour of classical pul- 
0_C sators was solved with the incorporation of turbulent convection into pulsation 
q hydrocodes. However, the reasons for the computed double-mode behaviour 
$_i ■ were never clearly identified. In our recent papers (Smolec & Moskalik 2008a, b) 
c/3 ' we showed that the double-mode behaviour results from the neglect of negative 
buoyancy effects in some of the hydrocodes. If these effects are taken into ac- 
count, no stable non-resonant double-mode behaviour can be found. In these 
proceedings we focus our attention on the role of negative buoyancy effects in 
k>" ' classical Cepheid models. 
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Introduction 

00 ' 

■ Since the early days of nonlinear pulsation computations, modeling double- 
, mode (DM) phenomenon was one of the major objectives. However, the search 
k>( ' f° r n on-resonant double-mode behaviour with radiative hydrocodes failed (see 
rS \ however, Kovacs & Buchler 1993). Typical modal selection observed, was, first 
overtone (10) pulsation at the hot side of the instability strip (IS), fundamental 
mode (F) pulsation at the red side, and either-or domain (F/10) in between. In- 
corporation of turbulent convection into pulsation hydrocodes led to stable and 
robust double-mode pulsation (Kollath et al. 1998, Feuchtinger 1998). Most 
of the so far published double-mode models were computed with the use of the 
Florida-Budapest hydrocode (eg. Kollath et al. 2002). This hydrocode adopts 
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time-dependent convection model based on the KuhfuB (1986) work. Kuh- 
fuB model was also adopted in our pulsation hydrocodes (Smolec & Moskalik 
2008a), however, we could not find a double-mode behaviour, despite exten- 
sive search for it (Smolec & Moskalik 2008b). We trace the difference in the 
computed modal selection to different treatment of negative buoyancy effects 
in both codes. In the Florida-Budapest code negative buoyancy is neglected, 
while it is present in our hydrocodes. Comparison of both treatments (per- 
formed with our hydrocodes) allow to understand the reasons for double-mode 
behaviour computed with hydrocodes that neglect negative buoyancy. 

Turbulent convection model 

In KuhfuB model of turbulent convection, equations of motion and energy con- 
servation are supplemented with additional, single equation for generation of 
turbulent energy, e t . Details of the model may be found e.g. in Smolec & 
Moskalik 2008a. Here we focus our attention on turbulent energy equation 
and its crucial terms (turbulent pressure and flux of turbulent kinetic energy 
neglected): 

% = S-D + B.V. (1) 

In the above equation, S is turbulent source function (or driving function), 
responsible for driving of turbulent eddies through the buoyant forces. Source 
function is proportional to superadiabatic gradient, Y = V — V a , and to the 

1 /2 

speed of convective elements, ~ e t ' . The D-term describes the decay of 

3 /2 

turbulent eddies, through the turbulent cascade, D ~ e t ' . E.V. describes 
energy transfer through eddy-viscous forces. It describes the interplay between 
turbulent motion and mean gas motion, being proportional to the speed of 
convective elements and to the spatial derivative of the scaled mean velocity 
field, U/R. It always damps the pulsations and contributes to the driving of 
turbulent energy 1 . 

The crucial difference between our and Florida-Budapest formulation is the 
treatment of turbulent source function in convectively stable (Y < 0) regions 
of the model. In our hydrocode we allow for negative values of source function, 
just as in the original KuhfuB model (S ~ Y). In Florida-Budapest approach 
source term is restricted to positive values only (S ~ Y + ), which is equivalent 
to the neglect of negative buoyancy. Following convention introduced in Smolec 
& Moskalik (2008a) we will denote the convective recipe and models ignoring 
negative buoyancy as PP models (Florida-Budapest approach), while convective 



1 Exact form of this term differ in different pulsation hydrocodes, sec Smolec & Moska- 
lik 2008a, however results presented here are insensitive to these differences. 



R. Smolec 



3 



recipe and models including negative buoyancy will be denoted by NN (our 
default formulation). 

Consequences of neglecting negative buoyancy 

Using our pulsation hydrocodes we have performed a detailed comparison of 
models, differing only in the treatment of the source function (NN vs. PP 
models). Crucial differences are observed for single-mode limit cycle (full am- 
plitude, monoperiodic oscillation) models. Amplitude of the models neglecting 
negative buoyancy (PP) is much lower than the amplitude of the models includ- 
ing negative buoyancy effects (NN). Lowering of amplitude in case of the PP 
models is connected with the eddy-viscous damping of pulsations in the deep 
convectively stable regions of the model. This damping is clearly visible in the 
nonlinear work integrals presented in Fig. 1. For PP model a significant eddy- 
viscous damping below zone 70 (marked with arrows in the Figure), not present 
in NN model is visible. These internal zones are convectively stable (Y < 0), 
however, in PP model significant turbulent energies are present in these zones, 
and hence eddy-viscous damping is possible. How these turbulent energies are 
built up? The bottom boundary of envelope convection zone, connected with 
hydrogen-helium ionization, is located roughly at zone 70 for both PP and NN 
models. Below this boundary, turbulent motions are effectively braked in NN 
model due to negative value of turbulent source function (eq. 1). Negative 
buoyancy slows down the turbulent motions very effectively and eddy-viscous 
damping is not possible in the inner parts of the model, due to negligible tur- 
bulent energies. Situation is different in case of PP models, in which turbulent 
source function is set equal to zero in convectively stable layers. Therefore, 
in these layers, turbulent energies are set by the balance between turbulent 
dissipation term (D term in eq. 1), which damps the turbulent motions, and 
eddy-viscous term, which drives the turbulent energies. Due to neglect of nega- 
tive buoyancy, turbulent motions cannot be braked effectively. To the contrary, 
they are driven at the cost of pulsation, through the eddy-viscous term. Below 
the envelope convection zone turbulent energies are as high as 10 9 — 10 10 erg/g 
- only three orders of magnitude smaller than in the center of the convection 
zone. Such high turbulent energies extend to more than six local pressure scale 
heights below the envelope convection zone, leading to significant eddy-viscous 
damping, visible in the left panel of Fig. 1. Deeper in the envelope, turbulent 
energies slowly decay, reflecting the vanishing amplitude of the pulsations. 

To check the effect of neglecting negative buoyancy on modal selection, we 
have computed several sequences of nonlinear Cepheid models using both PP 
and NN convective recipes. Details of mode selection analysis may be found in 
Smolec & Moskalik (2008b). Each model was initialized (kicked) with several 
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Figure 1: Nonlinear work integrals plotted versus the zone number. Surface at right. 



different initial conditions (mixtures of F and 10 linear velocity eigenvectors) 
and time evolution of the fundamental mode and first overtone amplitudes, A 
and Ax, was followed. Exemplary results are shown in Fig. 2. Computed trajec- 
tories are plotted with solid lines for PP model and dotted lines for NN model. 
Hydrodynamic computations coupled with amplitude equation analysis allow to 
find all stable pulsation states to which trajectories converge (attractors, solid 
squares in Fig. 2), and all unstable solutions, that repel the trajectories (open 
squares in Fig. 2). For single-mode solutions stability coefficients are computed. 
These are 7^0, which describes the stability of the fundamental mode limit cy- 
cle with respect to first overtone perturbation (switching rate toward 10) and 
70,1, which measures the stability of the first overtone limit cycle. Negative 
value of 7 means that the respective limit cycle is stable. If both coefficients, 
710 and 70,1 are simultaneously positive, both limit cycles are unstable, and 
double-mode pulsation is unavoidable. The run of stability coefficients across 
the instability strip, for sequence of Cepheid models, computed with both PP 
and NN convective recipes, is presented in Fig. 3. Arrow in this Figure marks 
the location of model, for which hydrodynamic integrations are presented in 
Fig. 2. 

Amplitude of the given mode is a main factor affecting its stability. Higher 
the amplitude of the mode, more able it is to saturate the pulsation instability 
alone, and hence, more stable its limit cycle is. For NN models amplitude of 
the fundamental mode is much higher than the amplitude of the first overtone 
(see Fig. 2) across significant part of the instability strip. Therefore, stability 
coefficient of the fundamental mode limit cycle, 710, becomes negative very 
close to the blue edge of the IS at temperature w 6290K (Fig. 3). Double-mode 
state is not possible, as first overtone becomes unstable (70,1 > 0) for much 
lower temperature (T « 6165K). In a temperature range in which both limit 
cycles are stable, either-or domain is observed, in which pulsation in either limit 
cycle is possible. 

In case of PP models, amplitudes of both modes are reduced as compared to 
NN models, but not by the same amount. As is well visible in Fig. 2, amplitude 
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Figure 2: Hydrodynamic integrations for a particular Cepheid model. Solid lines - 
model computed with PP convection, dotted lines model computed with NN convec- 
tion. 

of the fundamental mode is reduced much more than the amplitude of the 
first overtone. This effect is explained by higher amplitude of the fundamental 
mode (in comparison to 10) in convectively stable layers of PP model, in which 
significant turbulent energies are observed. As a result, eddy-viscous damping is 
stronger for the fundamental mode. Differential reduction of mode amplitudes 
is crucial in bringing up double-mode behaviour in PP model sequence. At the 
hot side of the IS, amplitude of the fundamental mode is lower, or comparable 
to the amplitude of the first overtone. Hence, fundamental mode limit cycle 
is unstable (71,0 > 0) in a wide temperature range. It becomes stable at 
much lower temperature (T » 6100K) in comparison to NN sequence. At 
this temperature first overtone is already unstable. Consequently, in relatively 
wide temperature range (~ 50K in Fig. 3) both limit cycles are unstable and 
double-mode state emerges. 

Summary 

Neglect of negative buoyancy has serious consequences for the computed Cepheid 
models. It leads to high turbulent energies in convectively stable layers, and con- 
sequently to strong eddy-viscous damping. This damping acts differentially on 
pulsation modes, which promotes the double-mode behaviour. If buoyant forces 
are taken into account, as they should be, no stable non-resonant double-mode 
behaviour can be found (Smolec & Moskalik 2008b). Therefore, the problem 
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Figure 3: The run of stability coefficients along a sequence of Cepheid models of 
constant mass/luminosity. Solid lines for PP models, dashed lines for NN models. 

of modeling F/10 double-mode behaviour in classical Cepheids remains open. 
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